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Abstract 



A sample of decays containing b- flavoured hadrons is tagged using leptons, and 
events having precise proper time measurements are selected. These events are used to 
study Bg oscillations. The flavour (b or b) at decay is determined from the lepton charge 
while the flavour at production is determined from jet charge or the charge of a second 
lepton, where available. The experiment was not able to resolve the oscillatory behaviour, 
and we deduce that the Bg oscillation frequency Amg > 5.2 ps~^ at the 95% confidence 
level. 
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1 Introduction 



The phenomenon of B — B mixing is now well established. In the case of the Bj system, 
the mass difference, Amd, between the two Bj] mass eigenstates has been measured rather 
precisely This mass difference gives the oscillation frequency between Bj^ and B^. Although 
these measurements can be used to gain information on the CKM matrix element Vtd, this 
is hampered by large theoretical uncertainties on both the meson decay constant, f^o, and 
the QCD bag model vacuum insertion parameter, B^o This difficulty may be overcome if 
Auig, the Bg oscillation frequency, is also measured. In this case, the CKM information can be 
extracted via the relation 

Amd m^o |FtdP /ro^bo 

where 777,30 and TTigo are the B^ and B^ masses, as the ratio of decay constants for and B^ 
mesons is much better known than the absolute values ^. Information on |Vtd| could then 
be extracted by inserting |Vts|, which is relatively well known 

The most restrictive of the published limits 0, |§, §|, 01 indicates that Arris > 9.6 ps^^ |^, 
while the previous best limit from OPAL gives Anis > 3.1 ps~^. 

A study of Bg oscillation is presented in this paper, using lepton charge to tag the flavour of 
the Bg at decay and a jet charge technique (Section 3), or opposite-jet leptons where available, to 
tag the Bg flavour at production. The analysis studies the oscillation as a function of the proper 
decay-time reconstructed using secondary vertices (Section 2), using a purity in semileptonic 
b-decays which is evaluated event-by-event (Section 4). Note that the analysis is sensitive to 
both Bg and Bj oscillation, but is optimised for the study of Bg oscillation; the parameters 
describing the B° system are input from previous measurements. 

The analysis technique is similar to that presented in previous papers [^, ^, but includes 
more data (from 1995), takes advantage of three-dimensional vertexing, and features a more 
sophisticated jet charge definition. Finally, it combines the single lepton and dilepton data. 



2 Event selection 



The analysis was performed on data collected by the OPAL detector at LEP in the vicinity of the 
Z° peak from 1991 to 1995. The OPAL detector has been described elsewhere 0|. Tracking of 
charged particles is performed by a central detector, consisting of a silicon microvertex detector, 
a vertex chamber, a jet chamber and z-chambers.Q The central detector is positioned inside a 
solenoid, which provides a uniform magnetic field of 0.435 T. The silicon microvertex detector 
consists of two layers of silicon strip detectors; the inner layer covers a polar angle range of 
I cos 6*1 < 0.83 and the outer layer covers | cos 6*1 < 0.77. This detector provided measurements 
of hits in the r-0 plane for data taken since 1991, with z coordinates also measured since 1993. 
The vertex chamber is a precision drift chamber which covers the range | cos^^| < 0.95. The 

^ The coordinate system is defined with positive z along the e~ beam direction, 9 and (f> being the polar and 
azimuthal angles. The origin is taken to be the centre of the detector and r is the distance to the z-axis. 
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jet chamber is a large-volume drift chamber, 4.0 m long and 3.7 m in diameter, providing 
both tracking and ionisation energy loss (dE/dx) information. The z-chambers measure the 
2;-coordinate of tracks as they leave the jet chamber in the range | cos^| < 0.72. The coil is 
surrounded by a time-of-flight counter array and a barrel lead-glass electromagnetic calorimeter 
with a presampler. Including the endcap electromagnetic calorimeters, the lead-glass blocks 
cover the range | cos6\ < 0.98. The magnet return yoke is instrumented with streamer tubes 
and serves as a hadron calorimeter. Outside the hadron calorimeter are muon chambers, which 
cover 93% of the full solid angle. 



Hadronic decays were selected using criteria described in a previous publication |10 
Only data where the silicon microvertex detector was fully functional were accepted, resulting 
in 2.9 million hadronic Z° decays with 3-dimensional silicon readout (3D), and 0.9 million 
with silicon readout only in the r-0 projection (2D). Tracks and electromagnetic clusters not 



associated to tracks were grouped into jets using a cone jet algorithm ITT 



The selection of electrons and muons as candidates for semileptonic decays of b hadrons 
followed the procedure described in the previous paper 0, the exception being that an improved 
purity tag at a later stage of the analysis allowed the cut on the lepton neural network akin to 
be loosened from 0.7 to 0.5. 



Monte Carlo events were generated using the Jetset 7.4 program with parameters tuned 



to OPAL data []I3|, and were processed by the detector simulation program |M. 



2.1 Proper time reconstruction 

An attempt to reconstruct the proper decay time of the parent b hadron was made for each jet 
containing a lepton with a minimum pt of 0.7 GeV/c relative to the jet axis. The algorithm 
used to reconstruct the decay length of the supposed parent b-hadron was different from that 
described previously [Q. In principle, each b jet should contain two vertices other than the 
primary vertex, corresponding to the decays of the b and c hadrons. These two vertices are 
referred to as the secondary and tertiary vertices, respectively. In the new algorithm, the 
positions of the secondary and tertiary vertices are allowed to vary independently in a maximum 
likelihood fit, where the likelihood is calculated for each vertex position by taking the product 
of likelihoods of all the tracks in the jet. 

Each track contributes a likelihood of {w/2) x (P^ + Pt) + {1 — w) x Pp, where w is the 
probability for the track to originate from the secondary vertex, as determined from the track 
momentum and angle relative to the jet direction. For the lepton candidate, w is fixed to 1. 
The quantities Pp, Pg and Pt are the probabilities for the track to be compatible with belonging 
to the primary, secondary or tertiary vertices respectively, based on the impact parameters 
relative to the assumed vertex positions both in r-0 and r-z. 

Such a fit has 6 free parameters: the B decay length, the D decay length and two angles 
each for the B direction and the D direction. In practice, the power of the fit to measure the 
D decay length was found to be poor. The resolution was therefore improved by imposing 
a constraint on the D decay length, L^: the likelihood was multiplied by exp(— L^i/L) where 
L = 0.093 cm was an average D decay length in the Monte Carlo, and Ld was constrained to be 
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positive. The effect of this constraint is to put the fitted D decay length to for most vertices, 
effectively reducing the number of free parameters to three. Gaussian constraints are imposed 
on the B direction angles, using the results of a B direction-finding algorithm similar to that 



described in [|15| , which weighted tracks according to their rapidity relative to the estimated B 
direction. Loose constraints were also imposed on the D direction angles. The B decay length 
was constrained^] to lie in the range -0.4 cm to 2.5 cm. 

For leptons coming from semileptonic B decays, according to Monte Carlo with 3D silicon 
information, about 4% of the jets considered were rejected as containing fewer than 3 tracks 
passing the quality requirements. In about 8% of cases the fit did not successfully converge, 
and these candidates were rejected. About 2% of the candidates were rejected with the fitted B 
decay length less than -0.4 cm. To ensure that the results of the fit were stable and unambiguous, 
candidates were rejected if the likelihood was improved when the decay length was increased by 
1 standard deviation. In addition, a scan of the likelihood was made as a function of the B decay 
length with the other parameters kept fixed. If the log likelihood was within 1.75 of the fitted 
result for any B decay length further than two standard deviations from the central value, the 
candidate was rejected. These two requirements rejected about 24% of the candidates. Thus, 
in total, about 62% of leptons from B decays were selected with a successful fit. For these 
leptons, the B decay length. Lb, and its uncertainty, ctl, were taken directly from the fit. 

Given a reconstructed secondary vertex, the B energy was determined in a similar way to 
that described previously @, |l6l. The energy of the jet containing the lepton was reconstructed 
using the Z*^ mass to constrain the event kinematics, and the estimated contribution from 
fragmentation particles was subtracted. The fragmentation particles were separated from the 
b-hadron decay products using momentum, angle and vertex information. The uncertainty of 
the B boost, a^-y, was determined from the estimated uncertainties on the charged and neutral 
fragmentation energies, which were determined from the estimated probabilities for each track 
or cluster to originate from a fragmentation particle. 

The selected vertices were split into two classes, A and B, according to the quality of the 
proper time reconstruction. The better quality secondary vertices (class A) were selected by 
requiring that 



• the mass of the reconstructed vertex, using association probabilities for each track as 
calculated in the likelihood, was larger than 0.5 GeV; 

• the per degree of freedom of the vertex configuration was less than 6; 

• the angle between the lepton direction and the vector joining the primary and secondary 
vertices was larger than 110 mrad; 

• the ratio of lepton pt (relative to the jet axis) and momentum p was less then 0.35; 

• the reconstructed jet energy was required to not exceed the constrained jet energy by 
more than 2 GeV. The reconstructed jet energy, computed as the sum of energies of all 
tracks and clusters in the jet, should normally be smaller than the constrained jet energy 
because of the missing neutrino from the semileptonic decay. 

^The constraint was implemented by subtracting a smooth penalty function from the log likelihood, which 
was in the quoted range and parabolic outside this range. 



6 
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Figure 1: The distribution of at for all selected vertices in the data (points), together with 
the Monte Carlo prediction (histogram). The vertical error bars are obscured by the size of the 
points. 



Class A contains 67% of the selected vertices for leptons from B decays in Monte Carlo; the 
remaining 33% constitute class B. 



The proper time, t, is given by 0: 



Lb 

Pi 



ttlb 



El 



(2) 



mi 



As in the previous analysis use is also made of the estimated uncertainty, cr^, on the proper 
time, calculated from the separately estimated uncertainties on the decay length, 0"^, and the 
boost factor, ap^\ 

where correlations between the uncertainties on Lb and /37 are neglected. This is unimportant 
because the shape of the t distribution is parametrised from Monte Carlo, where such corre- 
lations are included. The distribution of estimated cXf is shown for the data for all selected 
vertices in Figure |I] together with the Monte Carlo prediction. The slight discrepancy that is 
visible is not important, as the analysis uses the value of at estimated for each event. Plots 
indicating the proper time resolution for B^ vertices with estimated at < 0.3 ps are shown in 
Figure ^, together with curves from the fitted resolution functions, defined as in the previous 
paper Note that different resolution functions are used for class A and class B vertices. The 
plots are for the 3D silicon readout; the resolution functions were determined separately for 
events with 2D silicon readout, where the resolution is similar (though the efficiency is worse). 



We use the notation h — c= 1. 
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Figure 2: The proper time resolution for Monte Carlo B" vertices with estimated at < 0.3 ps 
separated into the two vertex classes, A and B, for two ranges of true proper time t'. The 
curves indicate the fitted resolution functions. 



8 



Events were accepted only if they contained at least one vertex with at < 0.3 ps. This 
proper time resolution is necessary to give sensitivity to Amg > 5 ps~^. A second identified 
lepton was considered only if separated by more than 60° from the first lepton (the leptons are 
ordered using p and pt). A second vertex in the event, containing a second identified lepton, 
was accepted if it had at < 0.6 ps. These requirements led to 43598 events selected in the 3D 
data, of which 5012 had a second identified lepton and 1788 had two valid vertices. For the 
2D data, 9452 events were selected, of which 1019 had a second identified lepton and 175 had 
two valid vertices. The distributions of reconstructed proper time are shown in Figure ^ for the 
single lepton and dilepton events together with the Monte Carlo prediction. Only one proper 
time per event, satisfying at < 0.3 ps, is included in the figure for the dilepton events. The 
proper time distributions are given separately for the 93-95 data and the 91-92 data (only 2D 
silicon information available). In both cases, the Monte Carlo gives a reasonable description of 
the data. Also indicated are the contributions from b — £ decays and from hadronic Z° decays 
to uu, dd, ss or cc. The remaining contributions from bb events, such as b — >• c — >• £ decays and 
background leptons, are not indicated, but are included in the total. 



3 Flavour tag 

In order to detect oscillation, one needs to determine the b-fiavour (b or b) of the 3° both at the 
time of production and decay. The fiavour at decay is inferred from the charge of the daughter 
lepton. The fiavour at production is obtained from the charge of a lepton, where available, in 
the hemisphere opposite to the first lepton. For events containing only one identified lepton, it 
is determined by a modified jet charge technique using information in both hemispheres. The 
rest of this section describes the modified jet charge technique. 

3.1 Lepton hemisphere 

The lepton hemisphere is defined by the axis of the jet that contains the lepton. For this 
hemisphere, an unweighted jet charge, as used in the previous analysis |Q, has the desirable 
property that the decay products of a neutral B meson contribute a net charge of zero, provided 
they are fully contained in the jet. The jet charge therefore depends only on the net charge of 
the fragmentation tracks, and hence depends on the produced b-fiavour. However, there is extra 
information that this method does not utilise. The leading charged meson from fragmentation 
should refiect the b production fiavour; also, it should have a different angular distribution and 
a different momentum spectrum than the other fragmentation products. For production, 
this leading fragmentation meson should have an enhanced probability to be a K meson. 

To exploit this supplementary information, a number of jet /hemisphere charge variables 
were defined: 

• the unweighted jet charge, J2i Qi, summing over all tracks in the jet; 

• the unweighted hemisphere charge, J2iQi, summing over all tracks in the lepton hemi- 
sphere; 

9 



T 1 1 1 1 1 1 1 1 1 1 1 1 1 

OPAL // 3D 




t [ ps ] t [ ps ] 

Figure 3: The proper time distributions for vertices with at < 0.3 ps for single and dilepton 
events with 3D and 2D sihcon information available (points) . The Monte Carlo predictions are 
superimposed (solid histograms), together with the contributions from b — > ^ decays (dashed) 
and from uu, dd, ss or cc events (dotted). 
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• Qhcm^ ^ weighted hemisphere charge, where a neural network is used to assign the weights 
(see below); 

• Qhem^ a variant on Qhem (see below). 

Each of these variables was multiplied by the lepton charge, Qi, to provide sensitivity to mixing, 
and fed into a neural network, trained to separate unmixed decays (where the flavour is the 
same at production and decay) from mixed 3° decays, with an output Qsame ('same' stands for 
the same hemisphere as the lepton). 

To determine the weights used for Q^ll^, a neural network was trained on hemispheres to 
distinguish fragmentation tracks whose charge is opposite to that of the produced b quark (as 
expected for the leading charged fragmention meson) from those that have the same charge. 
Four inputs were chosen per track: the rapidity relative to the jet, the track momentum divided 
by the jet energy and dE/dx weights for the vr and K hypotheses. B decay products were 
excluded from the network training. The neural network output, 6, gives a weak separation 
between right and wrong sign fragmentation tracks. A hemisphere charge was formed to exploit 
this: 

Qhem = EQ^x('^^-0.5) (4) 

i 

where the summation is over all tracks in the same hemisphere as the lepton. However, the 
output range with the best separation also includes a large background from B^ decay products. 
A second neural network, with output e, was therefore implemented to separate fragmentation 
tracks from B decay products using the track impact parameters and significances in both r-(j) 
and r-z projections. A variant on the above hemisphere charge was formed: 

Qhcm = EQ^x(^^-Mxe, (5) 

i 

where the summation is over all tracks in the same hemisphere as the lepton. 

This hemisphere charge is not optimal, because the important charge correlations between 
tracks (especially those from the Bg decay) have been neglected. In fact, when combined 
with the lepton charge, the unweighted hemisphere or jet charge gives a separation of mixed 
and unmixed B^ decays that is superior to that obtained with or Q^^. However, the 

correlation of Qhlm!''' with the unweighted hemisphere charge is not very strong, allowing the 
combination to achieve a superior performance. 

The neural network output for the combination, Qsame, is shown in Figure ^(a) for unmixed 
and mixed B^ decays. The power to separate mixed from unmixed B^ decays is approximately 
40% better than the unweighted hemisphere or the unweighted jet charge alone. The distribu- 
tion of Qsame for the sclectcd data events is compared to Monte Carlo in Figure ^a), where a 
reasonable agreement is seen. 

3.2 Opposite hemisphere 

The opposite hemisphere is defined by the axis of the highest energy jet in the event other 
than the lepton jet. For this hemisphere, a similar approach was taken, though the situation is 
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Figure 4: The distribution of (a) Qsame, (b) Qopp and (c) Qtot for selected Monte Carlo events 
where the lepton originates from a Bg decay. The components from unmixed and mixed Bg 
decays are indicated. 
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simpler in that it is not necessary to distinguish fragmentation particles from B decay products 
since both may carry useful charge information. A number of variables sensitive to the b 
production flavour were defined: 



i 

axis; 



Qiet = Z^i Q'i ( — ) where the sum is over all tracks in the highest energy jet in the 

\ ^bcam / 

hemisphere using /t = 1 and p\ is the longitudinal momentum of track i relative to the jet 



• Qjet using K = 0; 

• Qpvoh defined below, using a neural network output for each track in the hemisphere 
opposite the lepton; 

• vertex charge information - 2 parameters Qvtx and aq (see below). 

The charge variables were multiplied by the lepton charge and a neural network, trained to 
separate unmixed 3° decays from mixed 3° decays, was used to combine the 5 input variables, 
giving an output Qopp- 

For the Qprob variable, a neural network was trained to enhance tracks with the desired 
charge correlation with the lepton (opposite in the absence of mixing), using the same four 
inputs per track as for the track neural network used in the lepton hemisphere. A probability 
Q that the lepton parent was produced as a rather than a is formed for each track, i : 

Qi = Q^^ (C - 0.5) + 0.5 , (6) 

where Q is the output of this neural network for track i. The tracks are combined into an 
overall probability per hemisphere, assuming that the probabilities are independent : 

^.^. + ^.(l-^^.) ■ ^ ^ 

The power of Qprob to separate mixed from unmixed B decays, when multiplied by the lepton 
charge, is very similar to that of the conventional Q^^^ which was used in the previous paper 
The two quantities are strongly correlated, but some gain can still be achieved by combining 
them. Excluding the vertex charge information, the neural net combination of the 3 charge 
variables results in a hemisphere charge with a separation power 23% better than that of 
alone. 

Further improvement to the opposite hemisphere charge tag can be made by tagging charged 
B decays using vertex charge, as described in a recent paper Jl^. This gave two variables: Qvtx 
and ctq, which were available for about 40% of events (where vertices significantly separated 
from the primary vertex were found). This information, suitably transformed, was fed into the 
5 input neural network. The 5 inputs were still used when no vertex information was available, 
with the extra 2 inputs set to register the lack of information. A further improvement of about 
10% in separation power was obtained from the inclusion of the vertex charge information, and 
the network output, Qopp, is shown in Figure ^(b) for unmixed and mixed B^ decays. The 
distribution of Qopp for the selected data events, shown in Figure |^(b), is well described by the 
Monte Carlo prediction. 
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3.3 Combined flavour tag 



The Monte Carlo test samples indicate that the individual values of both Qsame and Qopp give 
the charge tag purity for Bg mesons fairly accurately. Since they are largely independent they 
may be combined simply to give a measure of the probability that the mixed: 

Q Qsame ^ Qopp X 2 ^ 

Qsame X Qopp + (1 " <5 same)(l Qopp) 

where the offset of —1 is introduced so that when Qtot = there is no information on whether 
the Bg mixed. The performance of Qtot is shown in Figure §(c), and has a separation power 
that is 40% better than Q2jet, used in the previous analysis |Q. The distribution of Qtot for the 
selected data events is shown in Figure ^, together with the Monte Carlo prediction. The good 
agreement observed serves as a check on the power of the charge tag. 



4 Tagging the b purity 

The selected lepton candidates contain significant contributions from sources other than the 
semileptonic decays of b hadrons, notably b — > c — > £ decays, background leptons either in 
Z — > bb events or in lighter quark events and c — > £ decays in Z — cc events . The proba- 
bility that a lepton candidate originated from a given source, ignoring lifetime information in 
the lepton jet, was evaluated using properties of the lepton jet and of the opposite hemisphere. 
These probabilities are used in the fit (described later) to scale the likelihood functions cor- 
responding to each source, where these likelihood functions describe the reconstructed proper 
time and charge correlation. 

For single lepton events, the value of akin, the lepton neural network output, was used 
together with a b-tag from the opposite hemisphere to produce a single variable that tags the 
purity of semileptonic decays of b hadrons. For dilepton events the purity was determined from 
the two values of Okin- 

The opposite hemisphere tag selects the highest energy jet in the thrust hemisphere opposite 
to that containing the lepton jet. In this jet, a neural network was used to separate tracks with 
good quality 3D impact parameter information coming from b jets from those coming from 
light quark jets (u, d or s). The inputs were the impact parameters and their significances 
in both projections, the distance of closest approach of the track to the jet axis and the 
for the track to intersect this axis. The neural network outputs for the individual tracks were 
multiplied together to form a jet probability, Pjso- Not all tracks with silicon information in 
the r-0 plane have reliable 3D impact parameter information, so a simpler variable was also 
constructed, based solely on the impact parameter significances in the r-cf) plane. In this case, 
Pj2D was constructed by multiplying the likelihoods for each track to be consistent with the 
primary vertex. The likelihood had a simple double- Gaussian form. A neural network combined 
PjSDi Pj2D, the numbers of tracks used for each of these quantities, the | cos^| of the jet axis, 
and Fox- Wolfram moments which also have distinguishing power between b jets and light-quark 
jets. The neural network was trained to separate b jets from u, d or s jets (charm events have 
properties that are intermediate). The neural network output, /^hemi, is shown in Figure ^ for 
bb events and events from lighter quarks. 
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Figure 5: The distribution of (a) Qsame, (b) Qopp and (c) Qtot for selected lepton events in the 
data (points) together with the Monte Carlo prediction (histogram). The dashed hne marks 
the value of the charge variable where there is no information on the b production flavour. 
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To form the event b-tag /3evt, a neural network was used to combine ctkin with /3hemi (though 
simply multiplying together the purities obtained from each tag would yield a similar perfor- 
mance). The distribution of /5evt is shown in Figure |^ for the data, together with the Monte 
Carlo prediction and the component from semileptonic decays of b hadrons. Also shown in the 
figure are the corresponding Pevt distributions with t > 1 ps. A comparison of these two figures 
indicates that the light quark background is reasonably modelled. 



5 Fit for Am 



s 



A maximum likelihood fit was constructed, similar to that described previously for single lepton 
events H and for dilepton events |]TB[. The following quantities were input per event for the 



three classes of events : 

1. single lepton events : t, at, Qtot and Pevt', 

2. dilepton events with one vertex : t, at, Q^l^ ■ Qf \ ctkin ^^"^ ^^"^^ 



kin' 



3. dilepton events with two vertices : t^'^\ at'\ t^'^\ af'\ Q^"* ■ Qf \ cx^L ^^"^ 



a 



(2) 
kin' 



where the superscripts are used to differentiate between the two leptons where appropriate. 
As mentioned before, the purities of the different sources were calculated event-by-event from 
the appropriate inputs, taking the distributions of these quantities from Monte Carlo for each 
source. 

The Monte Carlo predicts that bb events account for 81% of the single lepton events, 93.6% 
of the dilepton events with one vertex and 98.3% of the dilepton events with two vertices. For 
the selected vertices in bb events, 87.5% were predicted to come from b — > £ (or b — > c ^ 
b — >• r — > £, or b — J/V' — ^ ) Q decays and 9.2% from b — *• c — > £ (or b ^ J/^/; ^ £) decays. 
Of the b ^ £ decays, an estimated 10.5% involve Bg mesons. 

In this paper, the only parameter that is varied in the fit is the Bg oscillation amplitude. A, 
as defined in [|I^ and used in the previous paper [Q. The probability density for a produced 
Bg to decay as a Bg after time t' is: 

exp(-^) 1 + AcosAmst' 

' unmix ' Z 

r Z 

and to decay as a B^ is : 

_ exp(-^) 1 - AcosAmsf 

' mix ' n ' 

r 2 

where r is the B^ lifetime. At the true oscillation frequency, the fitted value of A should be 
consistent with 1. Far from the true frequency, the expectation value for A is 0. Therefore, 
values of Amg may be excluded when A is below 1 and inconsistent with 1. Since Amg is large, 

^ The decays b^c^^, b— >r^£ and half of the b — * J/t/i — > £^l~ decays are classed together with h ^ I 
because they have the same charge correlation between the lepton and the b quark. 
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Figure 7: The upper plot shows the distribution of (3evt for the data (points) together with the 
Monte Carlo prediction (solid line) . The lower plot shows the same thing but with a minimum 
reconstructed proper time of 1 ps. The components from the semileptonic decays of b hadrons 
(dashed line) and from uu, dd, ss or cc (dotted line) are also indicated. The latter component 
is almost invisible in the lower plot. 
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the fitted value of A is relatively insensitive to sources of systematic uncertainty that affect 
the mean Qtot or the overall like-sign fraction in the case of the dilepton events. Therefore a 
complicated fitting scheme, as used in the previous paper g], to reduce the impact of these 
systematic errors is not necessary here. This also implies that this analysis is not optimal for 
Arris values below about 2 ps~^. 

Unlike the previous analysis 0], the purity of the b production fiavour tag is taken to be 
independent of Arris or the fitted amplitude. The uncertainty on this tagging purity is assessed 
from the comparison of data and Monte Carlo shown in Figure ^. This agreement is also 
sensitive to the overall mixed fraction, x, averaged over all b hadrons - which is measured to 
be 0.118 ± 0.006 |]T| and the fraction of b — ^ c £ decays, taken to have a 15% uncertainty 

6 Results 

The performance of the fit was first tested using Monte Carlo simulated data. Figure ^ shows 
the fitted amplitude as a function of Arris for four samples generated with different true Arris 
values. These four Monte Carlo samples are statistically correlated, because the same Monte 
Carlo events are used to simulate different oscillation frequencies. Each represents 6.8 million 
hadronic Z decays (approximately double the statistics available in the data). The behaviour 
of these samples is consistent with expectation. At low values of true Arris, a clear peak is 
seen ai A = 1, while at higher frequencies the sensitivity is insufficient. In addition, tests were 
performed for true Arris values up to 15 ps~^ using additional Monte Carlo samples containing 
only Z — > bb events. The results of these tests, representing greater statistical power, were also 
consistent with expectation. 

For the data, separate results were obtained for the 91-92 data (2D silicon) (9452 events 
selected) and the 93-95 data (3D silicon) (43598 events selected). These results are shown in 
Figure ^, where the errors are statistical only. The data can also be split into single lepton 
(47109 events) and dilepton (6031 events) samples. These results are also shown in Figure |^. 
It can be seen that the single lepton results are more precise, except at very low frequency, 
where the dilepton events are sensitive to mixing on both sides of the event. No significant 
evidence for a signal at any particular frequency is seen in any of the four plots. The fact that 
the measured amplitude is consistent with at low frequencies, in each case, implies that no 
systematic problems are evident. 

6.1 Systematic errors 

Systematic errors arise from a number of sources. Their effect was determined by varying the 
appropriate parameter and obtaining a new set of amplitude results. This was done using the 
full data sample, with the exception of the resolution function uncertainties where Monte Carlo 
was also used, as described below. 

The sources of systematic error considered are given in Table |l]. In this table, fs and /baryon 
are the production fractions /(b B^) and /(b b baryon). The b — c — £ parameter 
represents a single scaling factor for the fraction of events due toh —>■ c —>■ i decays. The Z —>■ cc 
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Parameter 


Range 


/s 
/'baryon 

Z — cc 
Lepton background 


io.5r{:?% i 

10.1^1?% i 
±15% X nominal [g 
±20% X nominal § 
±20% X nominal 


Ts/Td 


1.564 ± 0.014 ps 
1.04 ±0.04 [|1J] 
0.99 ±0.05 il 
0.79 ±0.05 [| 


Amd 


0.464 ±0.018 ps-i [in 


Jet charge 
Resolution function 1 
Resolution function 2 


±0.01 in tagging purity at Qtot = 0.4 
worsen tracking resolution by 11% 
shift {xE)h by 0.02 



Table 1: Sources of systematic error. 



parameter represents a scale factor for the fraction of single lepton events due to Z — > cc decays. 
The effect is squared for the dilepton events. The lepton background number is a factor for 
the non-prompt background rate for each lepton. The quantity r+/rd represents the lifetime 
ratio for B"*" relative to B'' mesons, and similarly for Ts/t^ and t\^/tc^. These ratios also affect 
the composition of the sample, as the semileptonic branching ratios of the individual b hadrons 
are taken to be proportional to the lifetimes. As mentioned in the introduction, the analysis is 
sensitive to B^ as well as Bg oscillation, so a systematic uncertainty results from the value of 
Amd input. 

The systematic error due to the tagging purity of the b production flavour as determined 
from jet charge was mentioned above. The error was modelled by changing the relative separa- 
tion of the mixed and unmixed distributions of Qtot- The uncertainty on this offset was assessed 
using the comparison of the fraction of data with positive Qtot as a function of |Qtot| between 
data and Monte Carlo simulated events, taking into account uncertainties in the knowledge of 
X and the fraction of b — >• c — >■ £ decays. The table gives the size of the effect on the tagging 
purity at Qtot = 0.4 as an example. 

The resolution functions are determined from Monte Carlo (see Figure^, and so are affected 
by uncertainties in the simulation of the data. Two such effects are considered: firstly, the 
resolution of do and 0o was worsened for all tracks by 11%, where do is the impact parameter 
relative to the primary vertex in the r-0 plane, and (po is the (p angle at this point. This change 
represents the level of uncertainty in the tracking resolution. Secondly, the mean scaled energy 
of b hadrons, {xE)h, was lowered by 0.02. Such a change represents a shift of over 2a with 
respect to the measured value [0, but is inflated to include the effect of shape uncertainties in 
the b fragmentation function. Significant statistical fluctuations were observed in the estimated 
systematic errors due to the resolution functions. To reduce this effect, the systematic errors 
were estimated using a Monte Carlo sample together with the data, giving a total sample 
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size equivalent to 10.6 million hadronic Z° decays. Uncertainty due to charm fragmentation is 
expected to have a negligible effect, and was neglected. 

The individual systematic errors are given for the range of Arris values between and 
15 ps^^ in steps of 1 ps~^ in Table 0. The sign of the change in amplitude is indicated by ± or 
=F in each case, where ± indicates that A increases as the parameter in question increases. In 
the case of jet charge, the sign is defined relative to an increase in the charge tagging purity. 
For the resolution function uncertainty due to tracking resolution, the sign is defined relative 
to the change described above. The table also includes the overall fitted amplitude and its 
statistical error at these points. From the table it can be seen that the resolution function 
uncertainties dominate the systematic error at high frequencies, while the level of h c i 
decays is important at low frequencies. The combined amplitude results for the 91-92 and 93-95 
data together with the total errors are shown in Figure |10|. 



6.2 Exclusion regions 



To determine exclusion regions at 95% confidence level, we represent the measured value of 
A at a given value of Arris as a Gaussian distribution function G{A — /i,(T^), where /i is the 
central value and is the measurement error. Two alternative methods are then considered 
to determine whether the value of Arris is excluded: 



a) values are excluded where the probability of measuring an amplitude lower than that ob- 
served would be less than 5% were that value of Arris the correct one, i.e. 



/oo 
G{A- fi, a A)dA < 0.05 , (9) 



or 

b) the same definition, but limited to the positive region, i.e 

jrG{A-fi,aA)dA 



Jo^G{A-fi,aA)dA 



< 0.05 . (10) 



The first definition gives a true 95% confidence level, in the sense that there is a 5% probability 
to exclude the true value. However, it is not protected against setting hmits well beyond the 
experimental sensitivity. The second definition makes use of the fact that the predicted value 
of A lies between and 1, regardless of the value of Arris. It is automatically protected against 
setting limits beyond the sensitive range. For a true value of Arris well beyond the sensitive 
range, method (a) would exclude the true value in 5% of the experiments, while for method (b) 
this percentage would tend towards zero. 

For method (a) the excluded regions may be determined simply by plotting the curve A + 
1.645a A as a function of Arris (Figure 0). All regions where the curve lies below A = 1 may 
be excluded. This gives a lower limit Arus > 5.2 ps~^ at the 95% confidence level. The limit 
that would be obtained, were A measured to be at every value of Arris, is 7.0 ps~^. This is 
known as the sensitivity of the analysis. 
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Arus (ps ^) 





1 


2 


3 


4 


5 


6 


7 


A 


.02 


.05 


.16 


.31 


.26 


.25 


.49 


.75 


_stat 
'^A 


±.09 


±.16 


±.22 


±.27 


±.32 


±.39 


±.48 


±.57 


fs 


±.13 


±.16 


±.07 


±.02 


±.02 


±.01 


T.02 


T.05 


/baryon 


T.02 


T.Ol 


T.Ol 


T.Ol 


T.Ol 


T.Ol 


T.Ol 


T.Ol 




±.20 


±.24 


±.10 


±.04 


±.02 


±.00 


T.Ol 


T.Ol 


Z — > cc 


=F.02 


=F.04 


T.03 


T.03 


T.03 


T.03 


T.03 


T.02 


Lepton background 


±.05 


±.09 


±.10 


±.11 


±.11 


±.12 


±.13 


±.15 


(rb) 


T.OO 


T.OO 


T.OO 


T.OO 


T.OO 


T.OO 


T.OO 


T.OO 


r+/rd 


^.05 


T.Ol 


T.Ol 


T.Ol 


T.Ol 


T.Ol 


T.Ol 


T.Ol 




±.04 


±.03 


T.OO 


T.Ol 


T.Ol 


T.Ol 


T.Ol 


T.Ol 




^.01 


T.OO 


T.OO 


T.OO 


T.OO 


T.OO 


T.OO 


T.OO 


Amd 


±.05 


±.00 


±.00 


±.01 


±.01 


±.02 


±.02 


±.02 


Jet charge 


T.14 


T.18 


T.IO 


T.05 


T.04 


T.03 


T.06 


T.07 


Tracking resolution 


±.01 


±.02 


±.06 


±.06 


±.02 


±.00 


T.IO 


T.07 


b fragmentation 


±.00 


±.00 


±.05 


±.04 


±.08 


±.05 


T.Ol 


T.07 


syst 


±.30 


±.36 


±.20 


±.15 


±.15 


±.14 


±.18 


±.20 



Arus (ps ^) 


8 


9 


10 


11 


12 


13 


14 


15 


A 


.13 


-.07 


.48 


.66 


.43 


-.57 


-1.39 


-1.25 


_.stat 


±.69 


±.85 


±1.04 


±1.25 


±1.49 


±1.73 


±2.01 


±2.34 


fs 


±.04 


±.08 


±.03 


±.01 


±.07 


±.25 


±.39 


±.41 


/baryon 


T.Ol 


T.Ol 


T.Ol 


T.Ol 


T.OO 


T.OO 


T.OO 


T.Ol 




T.Ol 


T.Ol 


±.02 


±.02 


±.02 


±.05 


±.07 


±.12 


Z ^ cc 


T.03 


T.04 


T.04 


T.04 


T.06 


T.09 


T.ll 


T.13 


Lepton background 


±.15 


±.16 


±.21 


±.25 


±.29 


±.33 


±.37 


±.44 


(rb) 


T.OO 


T.OO 


T.OO 


T.OO 


T.OO 


T.OO 


T.OO 


±.01 


T+/Td 


T.Ol 


T.02 


T.03 


T.03 


T.03 


T.03 


T.04 


T.06 


rs/Td 


T.OO 


±.01 


±.04 


±.03 


±.03 


±.04 


±.06 


±.10 




T.OO 


T.Ol 


T.Ol 


T.Ol 


T.Ol 


T.Ol 


T.Ol 


T.02 


Amd 


±.02 


±.02 


±.02 


±.02 


±.02 


±.01 


±.01 


±.03 


Jet charge 


T.Ol 


±.02 


T.03 


T.05 


T.02 


±.05 


±.05 


T.Ol 


Tracking resolution 


T.09 


T.08 


T.Ol 


±.05 


±.25 


±0.97 


±1.43 


±1.59 


b fragmentation 


T.19 


T.31 


±.03 


±0.32 


±0.02 


T.14 


T.09 


TO.95 


syst 


±.26 


±.38 


±.23 


±.41 


±.40 


±1.07 


±1.54 


±1.96 



Table 2: Fit results for amplitude A, with the breakdown of systematic error contributions. 
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Figure 10: Combined amplitude results for the entire data. The errors shown include both 
statistical and systematic errors. Also indicated is the 'hmit' curve which is A + 1.645(7^. All 
regions where this curve lies below A — 1 are excluded in method (a). 
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For method (b), the confidence level must be calculated at each Arris point, according to 
the above formula. This gives the result Arris > 5.0 ps~^ at the 95% confidence level, very 
similar to the method (a) result. 



7 Conclusion 



Single lepton and dilepton events were used to study Bg oscillations with improved sensitivity 
with respect to previous OPAL papers 0, |^. The experiment was not able to resolve the 
oscillations, but can place a lower limit Arris > 5.2 ps~^ at the 95% confidence level. This result 
is consistent with previous publications [§, H, H, 0, and supersedes the previously published 
OPAL results |Q, The sensitivity of the analysis (the lower limit that would be expected 
were the true oscillation frequency very large) was found to be 7.0 ps~^, the second highest 
relative to previous publications 0] . 
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